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Abstract 



The behavior of the momentum distribution of top quarks in the threshold region is investigated. The 
quahtative behavior, in particular the dependence of the average momentum on the strong coupling 
O constant can be understood from analytical calculations for the Coulomb potential. Ambiguities in 
the relation between the excitation curve and the top mass are addressed. 
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CN ' Of particular interest for top quark studies in the threshold region is the dependence of the 
^ ' total and the differential cross section on the strong coupling constant. Some intuition and 
^ ; quahtative understanding can already be gained from the predictions based on a pure Coulomb 

■ potential. 

o • For a stable top quark of fixed mass the "effective threshold" can be associated with the 
o I location of the IS" resonance ^/sthr = '^^t + Eis with Eis = —E^yd = —a^nit/A which decreases 
^ ' with increasing a. The height of the resonance cross section is proportional to the square of 

■ the wave function at the origin and hence proportional to a^, as long as the the resonances are 
k>( ' reasonably well separated. In the limit of large F^, that is far larger than -E/j^^, the overlapping 

^ 15, 25* . . . resonances have to fill the gaps between the peaks. Since these gaps themselves 
■ - - ■ increase proportional to a^, one is left in the extreme case of large width with a cross section 
linear in a. Note that this corresponds to the behavior of the cross section close to but slightly 
above the threshold which is also proportional to a. 

For realistic top masses of about 150 GeV one thus observes a behavior of the peak cross 
section quite close to the first power in a. Since the location of the peak itself depends on a, 
only the analysis of the full shape allows to extract the relevant information. 

Once the threshold energy is determined experimentally, a and rrit are still strongly corre- 
lated and can hardly be deduced individually. 

Therefore in a next step also the momentum distribution of top quarks has to be exploited 
to obtain further information. The discussion is again particularly simple for the Coulomb 
potential V{r) = —a/r. The average momentum, in units of the Bohr momentum amt/2, can 
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be written in terms of a function /(e) which depends only on one variable e = E/Ef(yd if the 
energy E = ^/s — 2mt is measured in terms of the Rydberg energy. 

ip) = (1) 

For positive arguments the function / can be derived from obvious kinematical considerations. 

/(e) = v^ for e>0 (2) 

For the discrete negative arguments e„ = corresponding to the locations of the bound 

states the radial wave functions in momentum space are given in terms of the Gegenbauer 
polynomials 

HiPt = 1 . (4) 



with 



(2vr) 
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Using the explicit forms of 

Cl{z) = l, C\{z) = 2z, Cl{z)=Az^-l (5) 
one obtains through straightforward calculation 

= = /(-V9) = |^. (6) 

For arbitrary n one may exploit the following identity 

^, sin(n + l)0 



to derive the general result 
with the asymptotic behavior 



n^) (2n-l)(2n + l)7r 

/ (-4) - ^ . (9) 

This in accord with the result expected from classical mechanics: For the average momentum 
of a particle on a closed orbit in the Coulomb potential one derives 



, .n. ^ f<^\ - y ^ (1-e^cos^O" 

^ ^ - ^ '.i?/j,J 2vryo ^^1-ecosO^"-^ ' ^ ' 



Quantum mechanical orbits with angular momentum zero and high radial quantum numbers 
correspond to classical motions with excentricity e = 1. In this limiting case the classical ex- 
pectation value is easily evaluated and for n = 1/2 one finds agreement with the quantum me- 
chanical result. For small negative energies one therefore obtains the behavior /(e) = 2a/— e/vr. 
Significantly below threshold, however, the average momentum increases more rapidly with 
decreasing energy and between the 15" and the 25" state one observes an approximately linear 
dependence on the energy. 



2 



From these considerations the dependence of the average momentum on a (with E fixed) 
is easily understandable, in particular the seemingly surprising observation that well below 
threshold (p) decreases with increasing a. From (|I]) one derives for a shift in a (keeping the 
energy E fixed) the following shift in (p) 



'1-2^; 



a 



Above threshold as well as close to but below threshold / oc Y|e|. Hence ef / f = 1/2 and the 
average momentum remains unaffected. Significantly below threshold, however, e/'/ / ~ 1 and 
the factor in front of 5a/ a becomes negative. This explains why {p) decreases with increasing 
a. 

These results are illustrated in Fig.|l|. In Fig.^o we demonstrate that (p) as evaluated 
with the program for the Green function (solid line) coincides perfectly well with the values 
calculated from the analytical formula on resonance, indicated by the triangles. The prediction 
from classical mechanics, namely {p) oc \J\e\ is shown by the dotted line and agrees nicely for 
positive and negative energies. In Fig.|l|a as is increased from 0.20 to 0.24 and (p) changes in 
accord with the previous discussion. 

For definiteness we have chosen m = mj/2 = 60 GeV for the reduced mass and a = 4as/3 
with as varying between 0.20 and 0.24. To retain the stability of our numerical program 
Ft = 3 MeV has been kept non vanishing and a cut p < mt/2 has been introduced. The 
curves demonstrate the increase of (p) by about 10% for the corresponding increase in a. The 
triangles mark the locations of the resonances and the expectation values for the momentum 
as calculated from (^). 

The qualitative behavior remains unchanged for realistic QCD potentials. Predictions for 
different QCD potentials based on the techniques as described in p|, ^ are shown in Fig.||. 
These potentials correspond to different values of as{Mz)- Qualitatively the same behavior is 
observed as in Fig.|l|. In Fig.^a the top quark width has been set to an artificially small value 
of 0.03 GeV, in Fig.pD the realistic value of 0.3 GeV has been adopted. 

An important feature is evident from Fig.^: The momentum calculated for positive energy 
is nearly independent from and refiects merely the kinematic behavior, just as in the case of 
the Coulomb potential. This is characteristic for the choice of a potential where the large 
distance behavior is fixed by phenomenology and decoupled from the short distance value of 
as- 

At this point a brief discussion is in order on the relation between the top masses deter- 
mined on the basis of different potentials. The perturbative two loop QCD potential has been 
calculated in momentum space and is fixed unambiguously for sufficiently large Q^: 
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Figure 1: a) Average momentum as a function of E for different values of a. The markers 
show the results of the analytical calculation at l-S", 25", "iS energies, b) Comparison with the 
analytical result for discrete energies and with the square-root dependence close to threshold. 
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Figure 2: Energy dependence of (p), the average t quark momentum for ag = 0.13 (dotted) 
0.12 (dashed) and 0.11 (sohd) hue for rrit = 120 GeV. a) Tt = 0.03 GeV and b) Tt = 0.3 GeV. 
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Power law suppressed terms cannot be calculated in this approach. To evaluate the Fourier 
transform of this function in order to calculate the potential in coordinate space the small 
behavior has to be specified in an ad hoc manner. Different assumptions will lead to the same 
short distance behavior. The potentials will, however, differ with respect to their long distance 
behavior. In [§] it has been argued convincingly that the long distance tail is cut off by the 
large top width. However, an additive constant in coordinate space can be induced by the 
small momentum part of V{Q'^). This additional term would lead to a shift in the ti threshold, 
which in turn can be reabsorbed by a corresponding shift in m^. The different assumptions on 
the long distance behavior are reflected in differences between the predictions of |Q, ^ ^ for 
the precise location of the ti threshold for identical values of as and rrit and in differences in 
the as dependence of the momentum distributions for fixed rrit and energy (see also 0). All 
these differences can be attributed to the freedom in the additive constant discussed before. 
The same additive constant appears in hh spectroscopy, such that the mass difference between 
top and bottom is independent from these considerations. 
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